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Q_i' Abstract 

r~| , Using a toy model Lagrangian we investigate the formation of vor- 

tices in first order phase transitions. The evolution and interactions 
of vacuum bubbles are also studied using both analytical approxima- 
tions and a numerical simulation of scalar field dynamics. A long lived 
bubble wall bound state is discovered and its existence is justified by 
using a simplified potential for the bubble wall interaction. The condi- 
tions that need to be satisfied for vortex formation by bubble collisions 
are also studied with particular emphasis placed on geometrical con- 
siderations. These conditions are then implemented in a Monte Carlo 
simulation for the study of the probability of defect formation. It is 
shown that the probability of vortex formation by collision of rela- 
tivistically expanding bubbles gets reduced by about 10% due to the 
above mentioned geometric effects. 
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1 Introduction 



Symmetry breaking phase transitions in the early universe can give rise to 
topologically stable localized energy concentrations known as topological de- 
fects (for recent reviews see Refs. [0, 0, [3| ). These defects can be points 
(monopoles), lines (cosmic strings) or surfaces (domain walls) depending on 
the homotopy of the vacuum manifold of the broken phase 0. 

In general when a symmetry group G is spontaneously broken to a smaller 
group H, defects can form if the resulting vacuum manifold M = G/H has 
non-trivial homotopy. Cosmic strings (vortices in two space dimensions) form 
when the first homotopy group of M is non-trivial i.e. -K\(M) ^ 1. 

Consider for example the Lagrangian density describing the dynamics of 
a complex scalar field $ = $i + i$2 

£ = I^$W-V(|$|) (1) 

where V(|<&|) is minimized for |$| = a ^ (e.g. V($) = |(|$| 2 - cr 2 ) 2 ). The 
set of potential minima of V(|<&|) (<& = ae lOL ) has the topology of a circle S l . 
In a cosmological setup, according to the Kibble mechanism ||, there will be 
(by causality) field configurations that span the whole vacuum manifold as 
we go around a large circle in physical space (i.e. asymptotically $ — > ae td 
where 9 is the azimuthal angle in physical space). Such configurations will 
inevitably form || (but see Ref. @ for potential loopholes in the gauge 
case), with probability about | [||, when causally disconnected domains of 
the universe merge as the causal horizon expands. The asymptotic behavior 
$ — ► ae td implies by continuity of $ that there will be a point inside the 
large circle where $ = 0. This point (and its neighborhood) being outside 
of the vacuum manifold will be associated with topologically trapped energy 
density. This configuration is the topologically stable vortex ||. Extended 
to three dimensions this object becomes a line defect, the cosmic string. 

In the above simplified Lagrangian ([]]) no gauge fields are involved and 
the broken symmetry is a U(l) global symmetry leading to the formation of 
global vortices. The price to pay for considering a simple U(l) global rather 
than a U(l) gauge symmetry || is that the total energy of an isolated global 
vortex diverges logarithmically. This however is not a problem in systems 
where a physical cutoff scale is built in, like multivortex systems where the 
cutoff scale is the intervortex separation, or cosmological setups where the 
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cutoff is the horizon scale. For simplicity in what follows we will consider 
sytems with global rather than gauge symmetry breaking. 

The above picture of string formation when regions of the size of the 
horizon at the phase transition become causally connected, is characteristic of 
systems undergoing a second order phase transition ( []T0j| , see Jy] for a recent 



discussion). In systems undergoing first order phase transitions, vortices can 
form by the merging of expanding vacuum bubbles with scalar field phases 
such that the whole vacuum manifold is covered after the bubble collision and 
phase interpolation^. Consider for example an area AA of a two dimensional 
system undergoing a first order phase transition during a time interval AT. 
There are three basic conditions that need to be satisfied for a vortex to form 
within the AA during the time AT. 

1. The nucleation of at least three bubbles must take place during AT. 

2. The nucleated bubbles must have phases such that geodesic interpola- 
tion leads to complete coverage of the vacuum manifold S 1 . 

3. The initial geometric configuration of the three nucleated bubbles must 
be such that the collision of all three bubbles occurs before the phase 
interpolation process can spoil the previous condition. 

Thus, the probability for a vortex to form within AA during the time interval 
AT may be written as 

P tot (AA, AT) = B 3 P phase (B)S(B, v) (2) 

where B = T^^-AT is the probability for a true vacuum bubble of area 
A bub to form within AA during time AT, P p h aS e(B) is the probability that 
a geodesic interpolation of the bubble phases completely covers the vacuum 
manifold. The dependence on the probability B exists because P p h aS e{B) is 
larger for clustered defects and therefore it increases with B. The suppression 
factor S(B, v) depends on bubble formation probability and the bubble wall 
expansion velocity v, and is the probability that the initial configuration of 
bubbles will be such that the third condition is satisfied. 

The factor P p h aS e{B) has been calculated in previous studies and found 
to be between 0.25 and 0.42 depending on the number of clustered vortices. 

^ubcritical bubbles ]l2| can also play a minor role in vortex formation [ p"3| but those 
effects are ignored within our approximation 
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One of the main goals of this work is to find the suppression factor S(B, v) for 
relativistically expanding bubbles (v ~ 1). The case when bubbles expand 
in the presence of plasma |TJ| (v < 1) is significantly more complicated |IS 
and will be included in a separate publication I15J. 



In the next section we derive and solve numerically (and analytically in 
the thin bubble wall limit) the instanton equations for bubble formation. 
We use a simple scalar field potential with vacuum manifold S 1 describing a 
first order phase transition. The obtained scalar field configurations are then 
evolved by using a numerical simulation based on a second order accurate 
leapfrog algorithm [[H]]. We first focus on the bubble interactions and study 
their dependence on the phase difference between the two bubbles, both 
numerically and in the context of a simple analytic model. After justifying 
analytically the numerically observed existence of long lived states of repelling 
walls in the limit of large phase difference, we show that no metastable states 
exist in the model considered. Thus all colliding bubbles eventually merge 
and the existence of metastable embedded walls need not be a consideration 
for studying vortex formation. The formation of relativistic phase waves after 
bubble collisions is also seen in our simulations, confirming expectations from 
previous studies. 

In section 3, we focus on the probability of vortex formation during three 
bubble collisions. We first briefly discuss the factor P p h ase (B) and show how 
can its dependence on the vacuum decay rate T, which has been ignored 



in most previous studies, be used to explain previous numerical results [[13 



that seemed puzzling at the time of their derivation. We then focus on the 
geometric suppression factor S(B,v ~ 1) and first derive the geometrical 
conditions under which the vortex formation is suppressed. These conditions 
are then tested by using dynamical simulations of three bubble collisions, and 
a Monte Carlo simulation is constructed based on these conditions to obtain 
the dependence of the suppression factor S on the mean interbubble distance 
and therefore on the bubble formation probability B. Finally, in section 4 
we conclude and briefly discuss extensions of this work that are currently in 
progress 
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2 Bubble evolution and interactions 



Consider a complex scalar field $ = |$|e"*, in a (2+l)-dimensional space- 
time, whose dynamics is determined by the Lagrangian (|]) and a symmetry- 
breaking potential 



V = A 



l$l 



l$l 



(3) 



with e > (Fig. 1). 

The false vacuum ($ = 0) of the potential (3) decays via bubble nucle- 

])). The two 



ation to the true vacuuum (|$| = an = er|(3 + e + J (3 + e) 2 
dimensional field configurations of the bubbles nucleated during this first 
order phase transition can be obtained 
solving the Euclidean field equations 



T8 1 (see also [19j for a review) by 



d 2 $ 2d$_<9V(|$|) 
dp 2 p dp <9|$| 



(4) 



with p 2 = \x\ 2 + t 2 , t being euclidean time. The initial configuration of 
the field after tunneling has therefore an 0(3) symmetry. This symmetry of 
the initial bubbles will become 0(2+1) symmetry in Minkowski spacetime, 
where p 2 = \x\ 2 — t 2 . Solutions satisfy the boundary conditions 



as p 



oo ana — 
op 



as p — >• 



(5) 



Analytic solutions to (|J) can be found in the so-called thin wall approxima- 
tion, i.e. for e C 1, when the energy difference between the minima is much 
smaller than the height of the potential barrier. For the potential (3), the 
thin-wall solution is of the form PDl 



l$l 



a 
2 



1 — tanh 



a 



(p-Ro) 



(6) 



where R Q = 1/iyXae) is the bubble's initial radius, found by minimizing the 
total energy. 

In the general, or thick-wall case, solutions to (§) with the boundary 
conditions (^) can be obtained numerically using a relaxation technique |17] 
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An arbitrary constant phase is then asigned in the interior of each bubble 
and the resulting configurations are evolved by solving the dynamical field 
equations 

*-W=-f£ ( 7) 

where the dimensionless variables 

(j)=<&/a, x = xm, and t = tm (8) 

with m = V\a : were used. For this part of the code, a second order leapfrog 
scheme [T7| was implemented on a 400x400 lattice, using as a reference the 
algorithm of Ref. ||. Energy was conserved in all simulations to within 5% 
for the evolution timescales. 

Defining the bubble walls as p = R Q , we see that the 0(3+1) symmetry 
forces the bubbles to expand, approaching the speed of light on a timescale 
determined by the initial radius R D (|x& u b| 2 = R 2 a + 1 2 ). Because of Lorentz 
contraction, thick-wall bubbles rapidly become thin-walled ones. As a con- 
sequence, interactions between bubbles occur almost always in the thin-wall 
regime, but we have kept thick-wall bubbles as initial conditions for our sim- 
ulations for generality. For the figures presented, we have used e = 0.8 unless 
specified otherwise. 

It is evident from the simulations that the two-bubble interactions are 
strongly influenced by the difference between the field phases. Walls of bub- 
bles with phase differences Aa ~ tc tend to repell each other at short dis- 
tances, thus delaying the time of merging of the walls. As bubbles approach, 
this causes the walls to separate and come in contact again, producing an 
oscillating false- vacuum wall, as can be seen in Fig. 2. This domain wall 
eventually decays, its lifetime depending on the phase difference. The effect 
is clearly observed when Aa ~ 0.97T. For phase differences of Aa ~ tt, we 
have found lifetimes of 10, (in the dimensionless units of Eq. @). Points in 
the wall located in the line joining the bubble's center oscillate with a period 
T ~ 4 (e = 0.8). Notice that this oscillatory state occurs before the bubbles 
merge, and is not to be confused with the oscillations of the field's magni- 
tude that occur during bubble merging and after the phase interpolation has 
occured, as reported in |21[ and [22 1 . The origin of the effect in this case is 



the balance between the repulsive force originating from the phase difference 
between colliding bubbles and the attractive force of the scalar field poten- 



tial. On the other hand, in the oscillations discussed in Refs. pl| and p2| 
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only the attractive potential force is relevant since the phase difference on 
opposite sides of the oscillating false vacuum is zero. 

To better understand the numerically observed long lived oscillating states 
we consider the problem in the planar approximation and write the equations 
of motion of bubble walls located at x = ±£ (t) as 



dW 



(9) 



where a is the mass per unit area of the wall, and W the static interaction 
energy between the walls, per unit area. Writing the field as <fi = |0|e* a , we 
have (primes denoting derivatives with respect to the spatial coordinate) 



W 



dx 



A 2 



+ 



\ 2 (a') 2 



+ v{4>) 



(10) 



As a first approximation, consider V(4>), \4>\ and \<f>\' as constants in the region 
between the bubble walls and a' = Aa/2£. The force between the walls will 
be 

'* 2 (Aa) 2 



(ii) 



where and 0' are close but not excactly equal to zero. That is, the difference 
in phases produces a repulsive term dominating at short distances. This 
simple picture can be improved by considering the total field configuration 
as a sum of two bubble configurations, = |0|e ia , with |0| approximated by 
the sum of two thin-wall bubbles 



1 + tank 



x 



1 — tank 



'x + f 



(12) 



and a' = Aa/2£. The resulting potential energy W is plotted in Fig. 3 for 
different values of Aa and e, as a function of £. Reducing e has the effect 
of reducing the false vacuum potential energy, allowing the repulsive term to 
dominate at greater distances. The effect of increasing Aa is that of displac- 
ing the minimum of the potential to higher values of £. The approximation 
is expected to be valid only for £ ~ 1 i.e. before the bubbles begin to merge. 
Our simulations have shown that the gradient is actually locally time depen- 
dent and the initially binding potential of Fig. 3 eventually gets dominated 
by the attractive terms as the phase gradient decreases locally. 
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The existence of these oscillatory states leads naturally to the question of 
existence of "embbedded" meta-stable domain walls for this potential. Such 
domain walls would be static solutions to the field equations 



c/'+A|0|V = o (13) 

with boundary conditions — > as x 0; ^ <r?7 for x — > oo. A simple 
asymptotic analysis suffices to prove that there is no solution to ( pj|) satis- 
fying the boundary conditions. As x — > 0, we can always write |0| = C</ > x a , 
a = C a x b . However, inserting this in ( |TB"[ ) we obtain 

a = \ b=-l Cl = -\ (14) 

So the asymptotic form of the equations is not satisfied for a real coefficient 
C a , and therefore the solution with the required boundary conditions does 
not exist. Thus the only possible domain wall occurs for Aa = it, where the 
situation is analogous to that of a real field. Then Eqs. (13) become 

BV (15) 



with the boundary conditions = ±ar] at x — > ±oo, and = at x = 0. 
The solutions are the domain walls in this potential. For e nonzero, we can 
see that its existence is guaranteed by noting, as usual, that the problem is 
analogous to that of a particle moving in the potential — V(0), with the field 
representing the particle's position and the spatial coordinate interpreted 
as a time coordinate. The "particle" is energetically allowed to start from 
x = +ar] at t — — oo , and arrive at x — —at] at t — +oo. This is not the 
case when e = 0, and then only solitons that go from x = ±at] to x = are 
possible. But as we have seen, even in the case where this domain wall does 
exist, the asymptotic analysis above means that it is unstable under small 
variations of Aa (a static solution does not exist for Aa 7^ vr). As discussed 
above this was also verified by our simulations. 

After the two bubbles merge, the situation is the one already studied 



numerically in ]2l] and analytically in p^fl . A phase wave is generated in 
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the contact point, and expands inside the bubbles with the speed of light, 
interpolating between the original phases. The field magnitude oscillates in 
the region of contact, with amplitude that depends in the phase difference, 
this time being inversely proportional: a greater phase difference will result 
in a more energetic phase wave, that carries away the wall's energy more 
efficiently thus dissipating the vaccum oscillations rapidly. These oscillations 
however do not affect the phase interpolation: the phase waves are produced 
as soon as the interiors of the bubbles come into contact, and escape at light 
speed. 

3 Vortex Formation 

We turn now to discuss the probability of vortex formation in a three-bubble 
collision. As noted in the introduction, the probability of forming three bub- 
bles in an area AA within a period of time AT is B 3 . In order to form 
a vortex, the phase of these three bubbles must be distributed so that the 
interpolation following collision leads to complete coverage of the vacuum 
manifold (condition(2)). The probability for this to happen can be esti- 
mated as follows || : Consider a triangular lattice in a 2d physical space and 
asign a random phase to each point on the lattice. The question is What is 
the probability that a vortex lies within a given triangle? The mean phase 
difference of two neighbouring lattice points of a triangle is clearly 27r/4. 
In order that a vortex forms in that triangle, the whole vacuum manifold 
must be covered by interpolating the phases of the three lattice points and 
therefore the phase of the third lattice point of the triangle must lie on the 
opposite part of the phase circle. Thus, on the average, the phase of the third 
lattice point should be in a range of 2ir/4. This will happen with probability 
1/4. Thus, the probability for forming an isolated vortex in physical space 
is 1/4. In realistic cases however vortices do not form isolated but in vortex 
antivortex clusters especially when the vacuum decay rate is high. A very 
relevant question therefore is: What is the probability for forming a vortex in 
a vortex antivortex cluster?. The probability p + _ for forming an antivortex 
next to an already formed vortex is larger than the probability for forming 
an isolated vortex. The reason is that the mean phase difference between 
neighboring lattice points in a vortex surrounding triangle is not 27r/4 but 
27r/3 and therefore to form an antivortex next to a vortex we only need that 
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the fourth lattice point be in a phase range 27r/3. So p + _ = 1/3. Thus, the 
probability for having i antivortices around a vortex is 



Pi=i : \(p + -r(i-p + -y- i (i6) 




and the probability per defect in clusters of four or larger is p 3 /! = 0.43 which 
is significantly larger than the naively obtained result of 0.25. 

This effect of clustering explains the result found in Ref. | 13| , where a 
simulation of the vortex formation process is done by placing random-phase 
bubbles in a 2-dimensional lattice, and allowing them to evolve and collide. 
The author finds a probability of vortex and antivortex formation of 0.42 
(10 vortices and antivortices formed after the nucleation of 23 closely packed 
bubbles), instead of the expected 0.25, and atributes his result to unknown 
dynamical effects. However, if we notice that the chosen rate of nucleation 
produces a highly packed system, and calculate the probability using (Jl6|), 
we obtain 0.43, in agreement with the simulations. 

The dynamical delay of merging due to the phase repulsive potential, 
discussed in the previous section, does not affect vortex formation. The phase 
interpolation will occur eventually, and the crucial factor is not the time it 
takes, but the spatial range over wich the bubbles overlap, wich determines 
the spatial sectors where phases will interpolate. A phase repulsion (which in 
any case occurs for Aa extremely close to tt, an unlikely situation if the phases 
are required to be distributed in the vacuum manifold as in condition (1) of 
the Introduction), can only delay the interpolation event. In the formation 
of vortices, it is only important to know if (not when) the phases of each of 
the three bubble pairs will geodesically interpolate. 

In addition, for the case of bubbles with relativistic velocities (i.e., when 
the effects of the plasma are negligible), we have both the bubble walls and 
the phases propagating at v ~ 1. This means that the phase waves can never 
reach the bubble walls leading to a single new bubble with interpolated phase, 
before the third bubble has time to reach the two collided ones. 

Thus, in the absence of plasma, the important ingredient needed to find 
the probability of vortex formation is the initial geometric distribution and 
nucleation times of the colliding bubbles. In the case when the effects of 
plasma are ignored (v ~ 1) simulations show (Fig. 4) that we can picture the 
bubbles as circles, even when the merging has occurred, the phase wavefront 
continuing the circle formed by the walls. 
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For a vortex to form, then, it is easy to see that condition (3) of section 1 is 
satisfied iff the three expanding circles have at least one common intersection 
point. Thus, for example, the extreme case of three bubbles of the same size 
(i.e., nucleated at the same time), but with centers located along a straight 
line, will not form a vortex. If we allow for different nucleation times, bubbles 
will interact having different sizes, and then even non-aligned configurations 
will not lead to a vortex (Fig. 5). Notice that in contrast to what has been 
stated in previous studies |22| it is in general not possible to select a frame 
where all bubbles nucleate at the same time. The condition of simultaneity of 
three events is that there is a space-like planar surface that goes through these 
events. This is not always possible (consider for example the special case of 
three events on a straight line i.e. effectively in one spatial dimension). Even 
if it was possible to select a frame of simultaneity for three events in 2 + 1 
dimensions it would not be appropriate to do so in a Monte Carlo simulation 
as this artificial selection of a frame would introduce a bias in the measured 
probabilities. We therefore have used different times for the nucleation of 
each bubble. 

The equation of motion for the bubble walls is dictated by the initial 
configuration, as stated in the previous section, so that the radius of a bubble 
nucleated at time t = Ati will be r^t) 2 = R 2 Q + (t — Ati) 2 , the subindex i going 
from 1 to 3. Given the positions (xj, yi) and nucleation times of three bubbles, 
one has only to solve the system of three equations for the intersection point 
and time 

(x - Xi f + (y- y t ) 2 = R 2 + (t- At,) 2 (17) 

A vortex forms if the solution to the system is real, positive and finite. 

As a test for this condition, we have checked using the numerical simula- 
tion described in section 2 the formation of vortices and predicted formation 
times for several configurations. The geometrical condition of circle intersec- 
tion was found to be accurate in all cases, confirming the hypothesis that 
dynamical effects can be ignored in the case of relativistically expanding 
bubbles. 

Having tested the geometric model with dynamical simulations, a sup- 
pression factor for vortex formation can now be found using a Monte Carlo 
simulation. Random position and nucleation times were assigned for bubbles 
inside an area AA = (Ai? /2) 2 and configurations were subject to the re- 
quirement that bubbles do not form in an overlapping state. After rescaling 
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the considered area such that R — 1, we constructed 500 randomly chosen 
three-bubble geometrical configurations. The system ([H]) was solved in each 
case and the number of cases with no triple intersection was counted. We 
defined the ratio of the number of cases where triple intersection occured 
over the total number of cases as the suppression factor for vortex formation. 
This is the factor S(B,v ~ 1) of Eq. 2. A plot of S vs A ~ B is shown in 
Fig. 6. 

Clearly the geometric suppression factor is less important for high nucle- 
ation rate (small A) but it is not negligible for low nucleation rates. 

4 Conclusion 

We have studied the basic conditions that are needed for vortex formation by 
the merging of vacuum bubbles nucleated during first order phase transitions. 
There are three such conditions which include the existence of three colliding 
vacuum bubbles, the complete coverage of the vacuum manifold by geodesic 
interpolation of the bubble phases, and the existence of a triple collision point 
for the merging of the three circles that describe the relativistic expansion of 
the bubbles. The probability that each condition is satisfied was obtained and 
the result was compared with previous studies. In particular the existence of 
a triple collision point during the evolution of three colliding bubbles occurs 
with probability approximatelly 92% for low vacuum decay rate. Such a 
suppression of defect formation rate is not expected to modify in any major 
way cosmological models based on cosmic strings. 

We have considered the case of relativistically expanding bubbles and 
have therefore neglected the friction effects of plasma particles surrounding 
the expanding vacuum bubbles. Such particles, being massless in the false 
vacuum but massive in the true vacuum (inside the bubble), are expected to 
scatter on the bubble walls and decelerate them to non-relativistic velocities. 
In models where the effects of plasma are important enough to lead to slowly 
expanding bubbles, our analysis gives only an upper bound to the geometric 
suppression factor S(B,v) which is expected to rapidly drop as the bubble 
wall velocity v decreases. Indeed, for small bubble wall velocity v, the in- 
terpolating phase front propagating always relativistically inside the bubbles 
after the first collision, is more efficient in equilibrating the phases to the 
interpolated value before the third bubble reaches the bubbles that collided 
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first. Thus the formation of the vortex can be avoided much more efficiently. 
The study of the dependence of S on v in the presence of plasma requires the 
detailed numerical simulation of the effects of the plasma, and the construc- 
tion of a generalized geometrical model and Monte Carlo simulation based 



on these effects. That work is currently in progress [15 
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